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Abstract 

We investigate the existence and stability of solitons in an optical waveguide equipped with 
a Bragg grating (BG) in which nonlinearity contains both cubic and quintic terms. The model 
has straightforward realizations in both temporal and spatial domains, the latter being most 
realistic. Two different families of zero-velocity solitons, which are separated by a border at 
which solitons do not exist, are found in an exact analytical form. One family may be regarded 
as a generalization of the usual BG solitons supported by the cubic nonlinearity, while the 
other family, dominated by the quintic nonlinearity, includes novel "two-tier" solitons with a 
sharp (but nonsingular) peak. These soliton families also differ in the parities of their real and 
imaginary parts. A stability region is identified within each family by means of direct numerical 
simulations. The addition of the quintic term to the model makes the solitons very robust: 
simulating evolution of a strongly deformed pulse, we find that a larger part of its energy is 
retained in the process of its evolution into a soliton shape, only a small share of the energy 
being lost into radiation, which is opposite to what occurs in the usual BG model with cubic 
nonlinearity. 

PACS: 42.81. Dp, 42.65.Tg, 42.81. Qb 
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INTRODUCTION 



It is commonly known that solitons in various physical media are supported by balance 
between nonlinearity and dispersion. In particular, in the case of the nonlinear Schrodinger 
(NLS) equation, the cubic self-focusing nonlinearity, which is most typical in dielectric optical 
media Q, should be balanced by diffraction or anomalous chromatic dispersion, in order to 
provide for the existence of spatial or temporal solitons, respectively. If the relative sign of the 
diffraction/dispersion and nonlinearity is opposite, solitons (bright ones) cannot exist. 

The above arguments apply to intrinsic diffraction and dispersion in any optical material. In 
some media, however, the intrinsic chromatic dispersion is too weak to support solitons in the 
temporal or spatiotemporal domain; then, much stronger artificial dispersion can be induced 
by means of a Bragg grating (BG) written on the sample. Notable examples are temporal ||] 
and spatiotemporal |3j solitons that were created in second-harmonic- generating crystals. In this 
case, the size of available samples is a few cm, hence the corresponding soliton's dispersion length 
must be very small, < 1 cm, which cannot be provided for by the material's intrinsic dispersion, 
even if the soliton is very narrow, with the temporal width ~ 100 fs. As was experimentally 
demonstrated in the above-mentioned works [|| and ||] (and proposed in the theoretical work 
Q), the necessary strong dispersion can be generated by BG. A similar situation takes place 
in BG-carrying silica fibers, where very strong dispersion induced by BG makes it possible to 
generate solitons (supported by the usual Kerr nonlinearity) in a very short piece of fiber ||. 

The presence of BG gives rise to a dispersion relation between the frequency uj and prop- 
agation constant (wavenumber) k, to 2 = loq + k 2 , with the spectral gap < \ui\ < ujq, inside 
which families of gap solitons may exist Q. It can be readily seen that the expansion of the 
dispersion relation near k = yields lu^luo = ± k 2 / (2loq) + ..., i.e., it contains two branches with 
opposite signs of the dispersion. This simple argument implies that, unlike media dominated 
by the intrinsic material dispersion, there should be no essential difference between self-focusing 
(SF) and self-defocusing (SDF) nonlinearity when the dispersion is induced by BG. Note that 
a similar situation takes place for spatial solitons in a nonlinear planar waveguide, if strong 
artificial diffraction is induced by BG in the form of a system of parallel scores on a surface of 
the waveguide: in that case, the same dispersion relation is valid, with to and k realized as the 
longitudinal and transverse wavenumbers (see, e.g., Ref. 0). 

The lack of principal difference between the SF and SDF nonlinearities in BG systems suggests 
a possibility to introduce a model in which the dispersion (or diffraction) is induced by BG, 
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while the nonlinearity may change its sign with the increase of the input power. Then, one 
may expect that the model may give rise to two qualitatively different families of gap solitons: 
low-frequency ones, in which the SF nonlinearity is balanced by the dispersion branch with a 
sign corresponding to anomalous dispersion, and high-power solitons, supported by the balance 
between SDF nonlinearity and the normal branch of the dispersion. The simplest model of 
this type may be based on the cubic-quintic (CQ) nonlinearity, that has recently attracted 
considerable attention, chiefly in the multidimensional case, as the combination of the SF cubic 
and SDF quintic terms prevents collapse and makes it possible to anticipate the existence of 
stable solitons ||. It has been shown that CQ nonlinearity correctly describes the dielectric 
response of the PTS crystal Q. 

In this work, we introduce the BG-CQ model and study zero-velocity solitons in it. The 
model is meaningful in both temporal and spatial domains, the latter interpretation being most 
realistic. We find the solutions in an exact analytical form, and their stability is tested by direct 
simulations. We find, in accordance with the arguments given above, that the model gives rise 
to two disjoint families of stable solitons, one of which may be regarded as a generalization of 
the ordinary BG solitons in the model with the cubic nonlinearity [10], while the other family 
consists of essentially novel solitons, in which the quintic nonlinearity is dominating. 



THE MODEL AND EXACT SOLITON SOLUTIONS 



In the temporal domain, the model is based on a system of normalized coupled equations 
for amplitudes u(x, t) and v(x, t) of the two counterpropagating waves linearly coupled by the 
resonant reflection on BG: 

iut + iu x + [(l/2)\u\ 2 + \v\ 2 ] u-v [(1/4)M 4 + (3/2)M 2 M 2 + (3/4)|i>| 4 ] u + v = 0, (1) 
iiH - iv x + [(1/2)M 2 + \u\ 2 } u-v [(l/4)|w| 4 + (3/2)|v| 2 |n| 2 + (3/4)|u| 4 ] v + u = 0, (2) 

where the BG-induced coupling coefficient is set to be 1, and v > is a real parameter to 
control the strength of the quintic nonlinearity. Equations (|l]) and (||) assume the usual ratio 
1 : 2 of the coefficients in front of the self-phase and cross-phase modulation terms in the cubic 
part of the equation |IJ, and the ratio 1:6:3 for the quintic part was derived in Ref. [pd]]. 
The same equations apply to the spatial domain, i.e., a planar waveguide with BG in the form 
of a system of parallel scores, x and t being interpreted as the transverse and longitudinal 
coordinates. In fact, this realization of the model is most realistic, in view of the possibility to 
use the above-mentioned PTS crystal that gives rise to the CQ nonlinearity. 
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In this work, we confine ourselves to the study of zero-velocity solitons, which are sought for 

as 

u(x, t) = A(x) exp(i</>(x) — iujt), v(x, t) = B(x) exp(iip(x) — iujt), (3) 

where to belongs to the above-mentioned gap, 

- 1< w < +1, (4) 

and the real amplitudes A, B and phases (f), ijj satisfy equations obtained by the substitution of 
@ into Eqs. ([!]) and (||). After straightforward manipulations, two simple relations follow from 
those equations, 

±(A 2 -B 2 ) =0,-^(^ + ^=0, (5) 
ax ax 

For the soliton solutions, the fields A and B must vanish at x = ±oo, hence Eq. (||) yields 
B{x) = A(x). A constant value of (ft + tp can be set equal to zero by means of an obvious phase 
shift, hence we also have ift{x) = —<j)(x). Obviously, these relations between the amplitudes and 
phases imply that the two fields are subject to the relation u* = v. 

The remaining equations for the single amplitude A and single phase 4> take the following 
form: 

^ = 2asin(20), a = A 2 , (6) 
ax 

d(f> 3 5 2 

— = u) + -a--va +cos(20), (7) 

A quotient of Eqs. ([?]) and (^) yields an equation for <f> regarded as a function as a, 

dc c 3 5v ,„ ,. , oS 

— + - = -- + — a, c = w + cos(20), (8) 
da a 2 2 

which can be immediately solved to yield c(a) = co/a — (3/4) a + (5^/6) a 2 , where Co is an 
arbitrary constant. Only the solution with Co = is not singular, hence we finally have 

cos(20) = -uj - (3/4)a + (5^/6)a 2 . (9) 

The next step is to eliminate sin (ft from Eq. @ by means of the expression @ , arriving at an 
equation 



da / / 3 hu r . , 

= aJi-(„ + i .- 1 .«f] . do 



An integral of Eq. (|T^) can be written in an implicit form, 



, , r° db 

x(a)= I (11) 

26Jl-(u, + f&-§^) 2 



which describes soliton solutions in the interval < x < oo (for x < 0, the solution can be 
obtained trivially, as the function a(x) is even). Indeed, letting a — > in Eq. (|ll]), one has 
x — > oo, in compliance with the boundary condition that a(x) must vanish at \x\ —> oo. On the 



other hand, it also follows from Eq. (11) that a(x) attains its maximum value ao (which is the 
soliton's peak power) at x = 0. 

The integral in QTT| ) can be expressed in terms of incomplete elliptic integrals, but this formal 
representation is useless. It is more essential to find the soliton's peak power ao in terms of v 
and oj. Because da/dx = at soliton's center, ao must be a root of the right hand side of Eq. 
dloD, i.e., oj + (3/4)a — (5v/6)a 2 = ±1. This pair of quadratic equations have four roots; after a 
simple consideration, one can find that only two of them may be relevant, namely 



Oo)i 



_3_ 

5^ 



9 10i/ 
16 ~ IT 



(1-oj) 



(12) 



(«o) 2 = Vy 



3 9 10v„ . 



(13) 



Note that, in the limit u — ► 0, the root (|i~2| ) yields the peak power (4/3)(l — oj) of the Bragg- 
grating soliton in the usual model with the cubic nonlinearity |p.0| , while the root di~3| ) goes to 
infinity in the same limit. 



Further straightforward consideration of the structure of the exact solution (11) leads to a 



conclusion that, if the root (12) exists as a physical (real) one, i.e., in the region 



1-oj < 27/(160i/), 



(14) 



this root determines the amplitude of the soliton, while the root ( p~3|) is then irrelevant. Exactly 
at the border of the region (|i~4|), 1 — oj = 27/ (160^), no soliton exists, and in the region 



1-oj > 27/(160^) 



(15) 



the root (|l^) does not exist. Nevertheless, Eq. (|i~l]) yields a soliton in the region (|l5|) too, but 
it takes a principally different shape, with the amplitude given by the expression (|T^). 

Note that the solitons may exist only in the spectral gap (Q), or0<l — oj < 2, therefore 
the region ( [Tq ) is actually present only if v > 27/320 ~ 0.0844. Thus, in this case there exist 
two different families of the exact soliton solutions separated by the border 1 — oj = 27/ (160^). 
Typical examples of solitons of the two types are displayed in Fig. 1. In the region (|i~4|), the 
solitons may be regarded as obtained by a smooth deformation of the usual gap solitons known in 
the model with the cubic nonlinearity (T^j , while novel solitons existing in the region (|l~5| ) , where 



the quintic nonlinearity plays a dominant role, may have a characteristic "two-tier" structure, 
with a sharp (but nevertheless nonsingular) peak, as it is seen in Fig. 1. Another qualitative 
difference between the solutions of the two types pertains to their phase structure and parity: if 



the soliton's amplitude takes the value (12), it follows from Eq. (||) that cos(20(x = 0)) = — 1, 
and the value ( |l3| ) of the amplitude corresponds to cos(2cf)(x = 0)) = +1. It is easy to understand 
(with regard to the above normalization setting the constant value of (j)(x) + ip(x) equal to zero) 
that a consequence of this is that, in the solutions of the first type, Reu(x) and Rew(x) are 
odd, andlmu(x) and Imv(x) are even functions of x. In the solutions of the second type, the 
parities of the real and imaginary parts of the solutions are opposite. 

STABILITY OF THE SOLITONS 

It is known that, even in the model with the cubic nonlinearity, the soliton stability is a 
difficult issue. For the first time, a possibility of instability of a part of the gap solitons was 



predicted in Ref. [12] on the basis of the variational approximation (VA). In that work, three 
(quasi) modes of internal oscillations of a soliton were identified and each of which could become 
unstable, depending on the soliton's frequency to: one mode generated weak nonoscillatory 
instability below a critical frequency, to < ujH\ and the others gave rise to essentially stronger 
oscillatory instability in an interval — 1 < oj < u4r , with oj^ < u>il^ . Later, the stability of 
the gap solitons in the cubic BG model was studied rigorously by means of numerical methods 
applied to equations for small perturbations linearized about the soliton |0|. As a result, it was 
found that a part of the solitons are unstable indeed. Comparison with the predictions of VA 
shows that, while the above-mentioned weak nonoscillatory instability appears to be an artifact 
of the approximation (which can probably be explained by a general theory considering possible 
false instabilities predicted by VA for solitons H]), the stronger oscillatory instability sets in 
indeed if the frequency oj is smaller than a critical value, which is very close to that predicted 
by VA. 

By means of systematic simulations of Eqs. (0) and @j we have tested the stability of both 



families of the solitons given by the exact implicit form (11) against small initial perturbations. 
The resultant stability diagram in the (v, oj) plane is displayed in Fig. 2. The most important 
conclusion suggested by the diagram is that as nu increases the stability region for the usual 
solitons (corresponding to solitons in region (|l4|)) shrinks but the opposite occurs for the solitons 
of the novel type. 

The simulations also demonstrate that far from the stability border the unstable solitons 
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decay into radiation. However, in the vicinity of the stability border, after shedding some 
radiation, they rearrange themselves into stable solitons. Thus, with regard to the possibility of 
radiative losses, the stable solitons play the role of attractors in the present conservative model, 
i.e., the stable solitons are really robust objects. A typical example is shown in Fig. 3, where 
an unstable soliton evolves into a stable one. 

The model presents essential differences against the one with pure cubic nonlinearity not only 
concerning the existence of solitons and their stability against small perturbations, but also in 
dynamics of strongly perturbed solitons. The simplest and most important example of a strong 
perturbation is sudden uniform multiplication of the soliton's profile by an amplification factor 
a essentially exceeding 1, which is a result of the action of an optical amplifier. In the usual 
models with the cubic nonlinearity, it is well known that the amplified pulse separates into a new 
soliton and considerable amount of radiation. However, in the present model it is feasible that 
the quintic term, which has the sign opposite to that in front of the cubic one, may help to shape 
the amplified pulse into a new soliton, thus reducing the share of energy emitted with radiation. 
Figure 4, which displays typical examples of the relaxation of the strongly amplified (by the 
factor a = 2) pulses in the cubic model proper and cubic-quintic one, demonstrates that this is 
indeed the case. In fact, the effect of the quintic term (with v = 1/2) is very significant: in the 
cubic model, the final soliton retains only 11.6% of the initial energy, while the energy-retention 
share in the cubic-quintic model is 92.4%. 



CONCLUSION 



We have introduced a model whose linear part includes two counterpropagating waves coupled 
by the resonant reflection on a Bragg grating, and the nonlinear part includes cubic and quintic 
terms with opposite signs. The model gives rise to two different families of solitons (provided 
that the coefficient in front of the quintic term exceeds a minimum value, v = 27/320), which 
are found in an exact analytical form. One family may be regarded as a generalization of 
the usual Bragg-grating solitons supported by the cubic nonlinearity, while the other family, 
in which the quintic nonlinearity is dominant, includes "two-tier" solitons with a sharp (but 
nonsingular) peak. Also, the soliton families differ in the parities of their real and imaginary 
parts. In the plane (to, v), the two soliton families are separated by a curve on which solitons do 
not exist. Stability regions have been identified within each family by means of direct numerical 
simulations. It was also shown that in the vicinity of the stability border the unstable solitons 
do not decay. Rather, after losing some energy in the form of radiation, they evolve into stable 
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ones. In the present model, almost all the energy of a strongly perturbed soliton is retained in 
the process of its evolution, a fairly small share being lost with radiation, which is opposite to 
what occurs in the usual cubic model. 
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Figure captions 



Fig. 1. Typical examples of solitons produced by the exact solution (|ll|) in the regions (14) 
and ( |T5| ) for v = 0.135, shown, respectively, by the dashed and solid curves. The two examples 
pertain to w = —0.23 and u = —0.27 (both are taken close to the border between the two 
regions, 1 — u) = 27/(160^)). The inset shows a blowup of the latter soliton near its center, to 
demonstrate that it has no singularity at x = 0. 

Fig. 2. Stability diagram for both types of the solitons. The exact border between the two 
families is the curve 1 — u = 27/ (160^), see the text. To the left of the border (the region of the 
usual solitons), both stable and unstable regions are explicitly marked. The novel solitons for 
which the quintic nonlinearity is dominant occupy the region to the right of the border, where 
the white area represents unstable solitons. 

Fig. 3. An example of the evolution of an unstable soliton (with v = 0.1, u; = —0.1) into a 
stable one. 

Fig. 4. Evolution of a soliton pulse instantaneously multiplied by the amplification factor 
a = 2: (a) the usual cubic model, v = 0; (b) the cubic-quintic model, v = 1/2. In both cases, 
the original stable stationary soliton pertains to w = 0.8. 




Fig. 1. Javid Atai and Boris A. Malomed 
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